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Abstract 

We investigate steady state entanglement in an open quantum system, specifically a single atom 
in a driven optical cavity with cavity loss and spontaneous emission. The system reaches a steady 
pure state when driven very weakly. Under these conditions, there is an optimal value for atom-field 
coupling to maximize entanglement, as larger coupling favors a loss port due to the cavity enhanced 
spontaneous emission. We address ways to implement measurements of entanglement witnesses and 
find that normalized cross-correlation functions are indicators of the entanglement in the system. 
The magnitude of the equal time intensity-field cross correlation between the transmitted field of 
the cavity and the fluorescence intensity is proportional to the concurrence for weak driving fields. 



PACS numbers: 



I. INTRODUCTION 



The study of entanglement has emerged as a central theme of quantum physics in recent 
years. It is driven both by fundamental questions and by the increasing interest in applica- 
tions that go beyond the limit of classical physics. Entanglement as a measurable quantity 
is a complicated subject, in particular when the systems have multiple components. Here we 
choose to study entanglement and its possible avenues of quantification in an open quantum 
system. This system, the canonical model of cavity QED [1], has a single atom coupled 
to the mode of an optical cavity with two reservoirs or avenues for extracting information: 
spontaneous emission and losses from the cavity. 

Two particles (or systems), A and B are said to be in an entangled state if the wave 
function of the complete system does not factorize, that is \AB) 7^ |A)|.B). One consequence 
of this form of the wavefunction is that a measurement on system A yields information 
about system B without any direct interaction with system B. For systems with the same 
dimension, in particular, a (pure) state is said to be maximally entangled if tracing over 
one of the two systems, say A, leaves the other one in a totally mixed state; this means 
that one can gain complete knowledge of system B by performing measurements on A only. 
An example that is of relevance to this work is the maximally entangled state of an atom 
and a field mode, \^>) = (l/v^) (\^,g) + |0, e)) with the first index denoting the number 
of photons in the field mode and the second (e = excited, g = ground) denoting the state 
of the atom. A measurement of the state of the atom immediately tells us the number of 
photons in the field mode; or a measurement of the photon number immediately tells us the 
state of the atom. 

The von Neumann entropy E = —trA^PA^ogiPA) of the reduced density matrix of system 
A, pa = tr b{pab) Q] quantifies the amount of entanglement in a given bipartite quantum sys- 
tem in a pure state. For mixed states, on the other hand, although it is easy enough to define 
what is meant by a totally unentangled state — namely, one in which it is possible to repre- 
sent the density operator as an incoherent superposition of factorizable states — quantifying 
the amount of entanglement in a partially entangled state is not, in general, simple. The 
natural generalization of the pure-state measure indicated above, known as the entanglement 
of formation, utilizes a decomposition of the quantum state p = £V Pj\ipj)(ipj\ = ^ P3P31 
and then defines E = min(^2 • PjEj) where Ej is the von Neumann entropy for the density 
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matrix pj = and the minimum is taken over all the possible decompositions, which 

is in general a very challenging task 

as 

. As a result of this, alternative measures have 
been proposed, such as the logarithmic negativity [4!]. It is also possible that some partic- 
ular measurement scheme may result in a most natural unraveling of the density operator, 
in the sense of the quantum trajectories approach [5| (especially for systems that are con- 
tinually monitored), and in that case it may be physically meaningful to focus only on the 
entanglement of the (conditionally pure) states obtained via that particular unraveling. 

One of the main purposes of this paper is to determine how much information about the 
atom-field entanglement in our canonical cavity QED system can be gleaned from the kinds of 
measurements represented by the traditional correlation functions of quantum optics. As we 
shall show below, we are actually able to avoid the difficulties for mixed-state entanglement 
because, in the limit we are interested in, our system is, to a good approximation, in a pure 
state, in spite of its being an open system interacting with two reservoirs. 



II. CAVITY QED SYSTEM 

Fig. [1] shows a two level atom in a driven optical cavity. We consider a single-ended 
cavity, with the intracavity field decaying via the output mirror at rate k. The two-level 
atom has a spontaneous emission rate to modes out the sides of the cavity denoted by 7, 
which is generally less than the free space Einstein A coefficient. The resonant coupling 
between the atom and the field mode is given by g = p, eg y/u/2heoV with \i eg the electric 
dipole matrix element, uj the transition frequency and V the volume of the cavity mode. The 
driving field is taken to be a large classical field e incident on the input mirror, with small 
transmission T in , so that the incident flux (in photon units) inside the cavity is proportional 
to T in e 2 

The quantum trajectory wave function that characterizes the system under a non- 
Hermitian Hamiltonian is: 

00 

mt)) = J2 (C g ,n{t)e~ iE3 ' nt M +C e , n {t)e-^ t \e 1 n)) (1) 
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H = hg (a T <j_ + aa + ) — ina'a — i—a + a^ + ihe{a} — a) (2) 
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FIG. 1: Single atom in a weakly driven optical cavity. Here g is the reversible coupling rate between 
the mode of the cavity and the atom, k is the decay rate of the field mode of the cavity, 7 is the 
spontaneous emission rate, e is the external drive (taken to be a classical field). 



with collapse operators 
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associated with photons exiting the output mirror and spontaneous emission out the side of 
the cavity. The indices e(g) indicate the atom in the excited (ground) state, while n is the 
number of photons in the mode. The energies are E en = E g ^ n+ i = hw(n + 1/2). We have 
the usual creation (a') and annihilation (a) operators for the field, and Pauli raising and 
lowering operators o± for the atom. 

In the weak driving limit, the system reaches a steady-state wave function: 



where the Ajj are known 



*> = \0g) + A 1>g \lg) + A .e\0e) + A 2 j2g) + A 1>e \le) 
. They are 
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The quantities p and q would be 1 for coupled harmonic oscillators. In cavity QED they 
differ from unity due to the non-harmonic, or saturable, nature of the atom. The squares of 
coefficients of single excitation A lt9 , A Q e give the rates of detection of single photons through 
the output mirror or in fluorescence (steady state), while the squares of the double excitation 
coefficients A l e , give the rates of detection of two photons either in coincidence (one 
through the mirror, and one in fluorescence) or both out of the mirror. The variables are 

a = WT^) (10) 

(3 = ^-a (11) 

7 

p = 1 - 2C[ (12) 

(1 + 2CQ f . 

1 ~ (i + 2C 1 -2Cr i ) 1 ' 

C x = £- (14) 

K'J 

C[ = C 1 7 2K . (15) 
1 (2k + 7 ) K J 



The one-excitation amplitudes Ai y9 and Aq )£ are proportional to the driving field e; the two 



excitation amplitudes A 2 g , and A\ e are proportional to the square of the driving field, e 2 . 



0]. The norm of this wave function is ||^)| = v/T-p 0(e 2 ); hence to lowest order in e, the 
coefficient of the vacuum should be (1 — (l/2)0(e 2 )). The term 0(e 2 ) makes no contribution 
to lowest nonzero order in e for the correlation functions or entanglement measures considered 
here. 

The entanglement of formation for this system is calculated from the density matrix after 
tracing over the field variables: 



Patom = Tr field \*)m (16) 

(17) 



1 + A\ + A\ A l>e A, g + A c 



A 1>e A ljg + A 0je Aj e + A\ e 
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The eigenvalues of this matrix are, to lowest nonvanishing order, 

Ai = (A hg A , e - A he ) 2 
= \A hg \ 2 \A , e \ 2 (q-l) 2 

A 2 = 1 - (A l>g Ao, e - A lye f 

where q is defined in Eq. ( TT3|) . and we have defined 

2<? 
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The entropy E = — Ai log 2 Ai — A 2 log 2 A 2 is then (again to lowest leading order) 

< n4 
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where we have taken the weak field limit, e being the smallest rate in the problem, so 
tj k <C 1. The approximation ( 1221) will hold provided {t/ k) 2 <C |£|. 

This entropy is the same as that obtained by using the density matrix for the field alone, 
traced over the atomic degrees of freedon. 

The co,_ n rst induced by Wooters for two ^ ™ a,so be us ed to 
characterize entanglement between two quantum systems of arbitrary dimension 
111 ]. The concurrence for our system is 



C = ^2(1 -Trpl tom ) 
= ^A(A 1>g A Qje -A 1 , e ) 2 

= 2 (I) 2 1(| (23) 
To see why ]£] oc \A\ >e — Ao^Ai >g \ may be a good indication of entanglement, consider what 
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happens if the wavefunction is a product state. We could write 



|*)p = \iPf) ® \4>a) 

= (D \0) + Dxll) + D 2 \2)) <g> (C^) + C e |e)) 

= A)C S |0<7) + DiC^lly) + D o C e \0e) + D 2 C g \2g) + D 1 C e \le) (24) 

For weak excitations, the coefficient of the ground state of the system is D C g = 1, or 
C g — Dq — 1. Then the product state is 

|*) P = |0flr> + D x |l(/) + C e |0e) + D 2 \2g) + D l C e \le) (25) 

Just knowing the one excitation amplitudes does not yield any information about entangle- 
ment, as it is possible to have A lj9 = D\ and A 0jC = C e . A 2g gives no information about 
entanglement, just nonclassical effects in the field, as it only involves field excitation. For 
weak fields D 2 is exactly A 2g . The entanglement shows up in the value of Ai >e \ if this value 
does not satisfy A li£ = D\C e = Ao ;e A lg , then it is not possible to write the state as a 
product state. 

In the presence of a non-zero vacuum contribution (as any real quantum state will have), 
one can learn nothing about entanglement simply by measurement of one-excitation ampli- 
tudes or probabilities. For example, the state |0, g) + a(\l, g) + |0, e)) is entangled, but only 
if one is certain that the probability amplitudes for higher excitation are truly zero. A state 
of the form |0, g) + a(|l, g) + |0, e)) + 0(e 2 ) cannot be said to be entangled without informa- 
tion on the relative size of the probability amplitude A\ fi . Measurement of one-excitation 
amplitudes conditioned by a previous measurement can yield information about entangle- 
ment. This can be accomplished by utilizing cross-correlation functions. A first important 
conclusion out of this study is that a measure of the zero time cross correlation between 
the atom and the field, as well as the mean transmitted and fluorescent intensities yields a 
measure of entanglement in the weak field limit. 

III. ENTANGLEMENT FOR WEAK EXCITATION 

Equation (f22|) of the previous section gives the amount of entanglement in the system as 
a function of the one and two excitation amplitudes. In terms of specific system parameters 
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FIG. 2: A plot of C scaled by (e/7) 2 as a function of k/j and 5/7 for weak excitation, 
the concurrence is: 



This section analyzes the sensitivity of the concurrence to the different parameters that 
appear in Eq. (|26|) . while trying to give physical reasons for their influence on the entangle- 
ment. Despite the fact that the rates of decay could be the same through the two reservoirs, 
spontaneous emission (7) reduces entanglement more than cavity loss (k). This is due to 
the fact that a 7 event (spontaneous emission) must come from the atom, while a k event 
(cavity transmission) could come from either the drive or a photon emitted by the atoms 
into the cavity mode. A spontaneous emission event unambiguously leaves the atom in the 
ground state, and the system wavefunction factorizes. 

Fig. [2] shows a remarkable result in the entanglement of the system as a function of the 
three rates in the problem. There is an optimal value for the coupling constant g given 
a set of dissipation rates k, 7. For many interesting cavity QED effects, stronger coupling 
is generally better, such as the enhancement of the spontaneous emission by a factor of 
1 + 2Ci = 1 + 2g 2 //t7 (this formula strictly holds only in the bad cavity limit k » g, 7). 




(26) 
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However, here increasing the coupling of the atom and field mode eventually decreases the 
amount of entanglement. To explain this it is instructive to recall that the concurrence 
C = \2af3(q — 1)|, where a is the mean cavity field, and j3 = —ga/ r y is the mean atomic 
dipole. As the coupling g increases, for a fixed weak driving field e, the intracavity field 
a = e/ (n + 2g 2 /'j) decreases. The intracavity field is the sum of the driving field in the cavity 
e/n, and the field radiated by the atom, (— 2Ci/(l + 2C\))e/ k, the minus sign resulting from 
the fact that the radiated field is tt out of phase with the driving field on resonance. We see 
that as g and C\ increase, the intracavity field decreases. This means that the steady-state 
wavefunction has a larger vacuum component, and consequently less entanglement. Another 
way to view this is that the cavity enhancement of the spontaneous emission rate means a 
larger loss rate for the system as the coupling increases, which is bad for entanglement. 

More formally, consider what happens if the two-excitation amplitudes in Eq. ([5]) are 
arbitrarily set to zero, which amounts to setting q = in Eq. (1261) . in which case the 
entanglement is only determined by the prefactor \af3\. The steady-state wave function 
becomes 

M„ = |O|7> + «(|l0>--|Oe>). (27) 

7 

There are two interesting limits on this Eq. ( 1271) for the parameter f = g/j- If / ^> 1, 
the steady state wavefunction is approximately \ip) sa = |0)(|g) — fot\e)) which is a product 
state. Also, if / -C 1, the steady state wavefunction is approximately \ip) ss = |<?)(|0) + a|l)) 
which again is a product state. To have entanglement between the atom and cavity mode, 
we must have the parameter / ~ 1, so as to prepare a steady state wavefunction of the form 
\i/j) ss = I Op) + ac(\lg) — |0e)) = \0g) +a\—), a mixture of the vacuum with a small entangled 
state component. 

The decrease of the prefactor \a(3\ is the dominant reason why the concurrence decreases 
with increasing g for large coupling. Close inspection of Fig. [2] also shows that there is an 
optimal cavity loss rate k for entanglement for a fixed g and 7. This is a result of reaching 
a maximum in the population of the states different from the vacuum (Eq. (jSJ)). Our results 



here are consistent with the numerical results of Nha and Carmichael j^]. 

When the system is driven off resonance, its response is typically characterized by trans- 
mission and fluorescent spectra {12I . \\\ . Although these are important probes of the system, 
they do not, in this limit, carry information about the entanglement, since they are derived 
from only the one-excitation amplitudes. 
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FIG. 3: Contour plot of C as a function of 5/7 and A/7 for K/7 = 0.5 



The concurrence as a function of the detuning of the driving laser shows that the steady 
state entanglement decreases typically by a factor of 1/A 3 for large detuning, where A = 
(u — uj{) with uj the resonant frequency of the atom and cavity, and 07 the frequency of the 



driving probe laser. But in t 
at the vacuum-Rabi peaks 



re case where g is larger than k and 7, the response is maximized 



14| . Figure [3] shows a contour plot of C for parameters in the 



regime of cavity QED where the two decay rates are similar: 2^/7 = 1.0. The concurrence 
increases with increasing g on resonance up to a saddle point, and then decreases. However 
the entanglement persists for detunings on the order of g, the approximate location of the 
vacuum-Rabi peaks in the spectra of the system. 

Detuning to a vacuum-Rabi peak (A = ±g), generates a steady state wave function of 
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the form 

\iP) ss = |0, g) + aTi(jg/i)\l, ±) + a 2 T 2 (g/ 1 )\2, ±>, (28) 

where |n, ±) = (l/\/2)(\n, g) ± |n — l,e)) is the n photon dressed atom-field state one is 
tuned near and Ti(g/j) and T 2 {g/'y) are functions that are maximal when g ~ 7. This is a 
state of mainly vacuum, plus a part that has entanglement between the atom and the cavity. 
It would seem that by continuing to tune to a vacuum-Rabi peak as g increases, it would be 
possible to maintain the entanglement, but Fig. [3] shows that this is not the case. Rather, as 
argued (for the on-resonance case) above, the crucial parameter for maximizing entanglement 
is / = (7/7 oc l/y/n sat , where n sat = j 2 /8g 2 is the saturation photon number. This is the 
dependence on the nonlinearity of the atomic system. Recall that, if these were two driven 
coupled harmonic oscillators, q — 1 and there would be no entanglement. A nonlinear 
interaction between the two harmonic oscillators would be needed to entangle them, as in 
the signal and idler modes in optical parametric oscillation. This nonlinear interaction would 
generate two-mode squeezing, which could be measured by homodyne detection of mode 
A(B) conditioned on detection of a photon in mode B(A), just as squeezing in one mode can 
be detected via conditioned homodyne detection of a mode based on a photodetection from 
that modeQ, H]- The nonlinearity of the two-level atom is needed to generate two-mode 
squeezing and entanglement between the atom and the cavity field. Even though the driving 
field is weak and the atom never nears saturation, there can only be entanglement with a 
linear atom-field coupling if the atom has a nonlinear response, as two-level atoms do. 

The concurrence shows its sensitivity to different parameters. Fig. H] shows a contour plot 
of C versus (7/7 and A/7 for a case where the cavity decay rate is larger than the spontaneous 
emission rate (k/7 = 10.0). The entanglement is largest near g/j — 4.0, before the vacuum- 
Rabi splitting of the spectrum, which does not occur in this case until (7/7 ~ 10.0, at 
which point the entanglement is already diminishing. The size of the maximum concurrence 
decreases by increasing ft/7 from 0.5 to 10.0 by a factor of about 30. 

IV. MEASUREMENTS OF ENTANGLEMENT WITH CORRELATION FUNC- 
TIONS 

The calculation of entanglement leads now to the question of how to implement mea- 
surements that give the full information in the case of this cavity QED system under weak 
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FIG. 4: Contour plot of C as a function of 5/7 and A/7 for /c/7 = 10 



excitation. The previous section shows that the concurrence is related to the rate of sin- 
gle photon counts out of the cavity or in fluorescence and to the rate of coincident counts 
from the cavity and fluorescence. These are the quantities associated in quantum optics 
with correlation functions, first introduced by Glauber [l?], Q, Q, Q. Generally these 
correlation functions involve comparing a field (intensity) of one mode with the field (inten- 



sity) 
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;he same mode at a later time (or different spatial location), with some exceptions 
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261 ] . However, entanglement in cavity QED has two components: atom 
and cavity mode. It is natural to look at cross correlations between the cavity mode and 
the fluorescent light that falls in the mode of the detector. 
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Consider a general cross-correlation function for two-modes of the electromagnetic field: 

G= (h{b\b)f 2 {a\a))/{h{b\b))(h{a\a)). (29) 

with fi and fa well behaved functions, in the sense of a convergent Taylor series on the 
Hilbert space of interest. If \ifj) is a product state, the correlation function G(a, b) factorizes 
and then is unity. If it is not a product state, then this will manifest itself in a non-unit 
value for the normalized cross-correlation functions. 

The simplest cross correlation function to consider is g^i^O)- This could be obtained 
by measuring the visibility of the fringe pattern formed by interfering the transmitted and 
fluorescent light. For the weakly driven cavity-QED system, this is 



(a) 
a(3 

= 1 (30) 

so to lowest order, there is no information in this correlation function about entanglement. 

To obtain information about entanglement the correlation function has to probe the 
two-excitation part of the state. A possibility to do this is the intensity cross correlation: 

(2) (o-+aW_) 



{a)a){o- + oJ) 



\A h 



|2 



l^lg^Oel 2 



= q 2 (31) 

This normalized correlation function is directly related to the coefficient of double exci- 
tations (See Eqs. (J5j), (JE]), (jTBl ). If q — 1 then fi'^(O) = 1 and there is no entanglement; 
so a non-unit value of q indicates entanglement. Using second-order intensity correlations 
las been proposed in the context of entangled coherent states by Stobihska and Wodkiewicz 



27] . 



(2) 

The cross correlation function 5^(0) contains information about the average photon 
number in coincidence with a measurement of the fluorescence relative to the average photon 
number in the absence of any interrogation of the fluorescence. ^^(O) — 1 = q 2 — 1 is an 
indicator of entanglement. 
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A way to measure g directly utilizes a field- intensity correlation function hg(r) [281 ]. 
that can be implemented as a homodyne measurement conditioned on the detection of a 
fluorescent photon, 



7 /■ , rt . (IfEt) 



^o(O) 



(If)(Et) 
((at + a)<j + a_) 
(at + a)(a + a_) 

= q (32) 

So hJ= Q (0) — 1 = q — 1 is also an indicator of entanglement in this system. What makes 
this measurement possible experimentally is the conditioning that selects only times when 
there is a fluctuation and the rest of the time (when the vacuum is present) no data is 



29]. 



collected |16|]. For one mode, the homodyned transmitted field conditioned by detection 
of a photon from that mode, is a measure of squeezing in that mode [28| . A homodyne 
measurement of the transmitted field conditioned by detection of a fluorescent photon is 
a measure of the two-mode squeezing, with the cavity field and atomic dipole as the two 
components. Generally, two-mode squeezing is an indicator of entanglement between the 
two modes. Gea Banacloche et al. explored this correlation function in a different regime of 
cavity QED and found it to be a witness of the dynamics of entanglement 

Non-classicality and entanglement are not necessarily simultaneously present. For exam- 
ple for two oscillators one could have = (1/ V2)(\A, B) + \B,A)), where A and B are 
coherent state amplitudes. In this state, there is entanglement, but each individual mode 
shows no non-classical behavior. Conversely, one can have non-classical behavior with no 
entanglement, say for example the atom in the ground state and the field in a squeezed 
coherent state. 

There is a particular form of the Schwartz inequality that must be satisfied for a classical 
field for the specific case of the system we are considering here: 

0©(O) - l) 2 < IG7&0) - 1)0^(0) - 1)1, (33) 

Here TT and FF denote zero delay intensity correlations for the transmitted and fluorescent 
fields respectively. In the one-atom limit, gp F (0) = 0, and (^(O) = q 2 p 2 , so this inequality 
becomes \q 2 — 1| 2 < \q 2 p 2 — 1] which depends on q, but also on the parameter p (Eq. f|T2l ). 
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which can be varied independently. There is no one-to-one relationship between Schwarz 
inequality violations and entanglement (by this measure) in this particular system. 

V. CONCLUSION 

We find that entanglement in weakly driven cavity QED is characterized by compari- 
son of two-excitation probability amplitudes to single excitation amplitudes, in particular 
the amplitude involving one excitation in each subsystem. It is necessary to have a small 
saturation photon number to enhance the nonlinear response which generates a larger en- 
tanglement. But this is true only to a point. We find the maximal entanglement for small 
k and when g/j is on the order of unity. This stems from the dual role of the coupling g. It 
couples energy into the atom, but due to cavity enhanced spontaneous emission, it can also 
channel energy out. 

Increasing 7 decreases the entanglement, and this can be explained in terms of the effect 
of the two decay processes on the system. If we detect a fluorescent photon we know it 
has come from the atom, and the atom is in the ground state. If we obtain a transmitted 
photon, it could have been emitted from the atom into the cavity mode, or just be a driving 
field photon that has passed through the cavity without interaction with the atom. It is the 
interference of these two indistinguishable processes that leads to nonclassical effects in the 
transmitted field. 

We have found a variety of cross-correlation functions that are indicators, or witnesses, of 
entanglement in this system. One can learn nothing about the entanglement by examining 
only one- or two- excitation amplitudes separately. In particular we find that a measurement 
of two-mode squeezing, or a homodyne measurement of the transmitted field conditioned on 
the detection of a fluorescence photon is directly proportional to the entanglement calculated 
via the reduced von Neumann entropy. Further work remains to generalize this approach to 
situations with higher drives, but the general approach of looking at entanglement together 
with the specific correlation function to measure gives physical insight into this problem. 

We would like to thank J. P. Clemens for fruitful discussions related to the topic of this 
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